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We present a new nonperturbative method to deal with the time-dependent quantum many-body 
problem, which is an extension of Wegner's flow equations to time-dependent Hamiltonians. The 
formalism provides a scaling procedure for the set of time-dependent interaction constants. We 
apply these ideas to a Kondo model with a ferromagnetic exchange coupling switched on over a 
' time scale r. We show that the asymptotic expectation value of the impurity spin interpolates 

continuously between its quenched and adiabatic value. 



> 

o 



- 1—1 

x 



PACS numbers: 72.15.Qm 



I. INTRODUCTION 



D ' The time-dependent quantum many-body problem has become of particular interest within recent years. New 
5^ . experimental investigations offer many different possibilities to study quantum systems in the presence of external 
driven degrees of freedom, such as cold atomic gases [l| , ultrafast spectroscopy on semiconductors @ , flux pattern 
formation in superconducting films @ or single molecule force spectroscopy [j]. Many of the available theoretical 
methods evaluate the quantum dynamics of such driven systems in a certain limit. If the time-scale of the external 
degrees of freedom is slow with respect to the intrinsic dynamics of the system, it is possible to perform calculations 
with high mathematical rigor by means of adiabatic perturbation theory Q. Many methods apply to time- 
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periodic perturbations. If the driving frequency lu exceeds the largest energy-scale of the unperturbed system, 
good method is an expansion in powers of g/ui, with g the coupling between the system and the perturbation, to 
' derive a time independent effective Hamiltonian Q. In the general time-periodic case one often follows the Floquet 
approach Q. The problem of calculating the quantum evolution of a certain initial state for an arbitrary switching 
procedure is often addressed in time-dependent Keldysh perturbation theory js}. However, such calculations are 
frequently limited to linear response, and expansions with respect to bare couplings can produce errors which grow in 
time. Furthermore, self-consistency schemes often yield a system of coupled non-linear integro differential equations, 
which are not convenient for numerical evaluation. Pure numerical techniques are the time-dependent density-matrix 
£N) . renormalization-group [10| , time-dependent density- functional theory jllfl , and exact diagonalization of the operator 
t-H ' H — idt in discretized form. 

I . Here we develop a new method to deal with the driven many-body problem, which is a similarity renormalization 
scheme for time-dependent Hamiltonians based on Wegner's flow equation method for the time-independent case 
[L?! [l3j . We apply our new approach to the isotropic ferromagnetic Kondo model 



H(t) = J2 £kc{ a c ka + J2 J(t) S ■ s kk , (1) 

k,a k,k' 

in a situation where the ferromagnetic exchange interaction J(t) < is switched on continuously from zero to a 
nonzero constant value Jo- Here S is the impurity spin-1/2 degree of freedom and Skk' = \ Y] n a c\ a <5 a pCk' ,p is the 
conduction band electron spin density. We choose this example for illustrating our approach because it is a many-body 
system which allows for a perturbative evaluation of the truncated RG equations, which posess a nontrivial RG flow. 

The second motivation comes from recent calculations that provided insights in the quantum dynamics of systems 
subjected to a sudden interaction quench \\W Il5j . These systems exhibit nonequilibrium dynamics in which the 
expectation values of certain observables are enhanced with respect to their equilibrium values. Specifically, we 
consider a time-dependent Hamiltonian 

H(t) =H + f(t) Htm (2) 

where the term H- lnt is switched on by some protocol fit), fit — > — oo) = 0, fit — » +oo) = 1. We study the time 
evolution of the system prepared in the ground state |Qo) of Hq in the infinite past, that is |Oo(^)) generated by H it) 
with \Cto(t — > — oo)) = Let O be any observable that commutes with Hq. We define its mismatch factor /i 

dcf (O)ncq ~ (O)o , s 

" (0) cq - (O) [6} 



Here (O)o — {Q,q\0\Qo) is the noninteracting expectation value, (0) cq = (fi|0|fi) is the interacting expectation value 
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with respect to the ground state of H = Ho + H lnt and 

(o) ncq = hm 1 [ T dt (n a (t)\o\n (t)) , (4) 

is the time-averaged nonequilbrium expectation value. 

For adiabatic switching on of i?j n t we have \i — 1 since | £^o (*)) asymptotically becomes the interacting ground 
state |fi). Remarkably, in the opposite limit of sudden switching f(t) = Q(t) of a weak perturbation H int one finds 
the universal value /j, = 2 in many quantum systems. In fact, a mismatch factor fj, = 2 can be proven rigorously 
for discrete weakly interacting systems [16| in second order perturbation theory. The same factor \i = 2 arises in 
many weak-coupling continuous systems like the quantum sine-Gordon model for 1 > 8ir [17| or in the ferromagnetic 
Kondo model with an infinitesimal positive magnetic field [H], [HI- In the latter model it describes the ratio of the 
large time asymptotic impurity spin normalized through its equilibrium value 

_ (S^(t->oo)) ncq -l/2 
fI - <S,) oq -l/2 (5) 

For an interaction quench in a Fermi gas in d > 1 dimensions, the mismatch factor of the correlation-induced 
corrections to the momentum distribution between nonequilibrium and equilibrium state turned out to be the key to 
understanding thermalization [l5l [l6j. However, in this model the universal value /i = 2 is only found when averaging 
over an intermediate well-defined time window in (j4|), which describes the prethermalized regime. 

Due to the important role of the mismatch factor fi in characterizing the nonequilibrium dynamics, it is natural to 
investigate its behavior when the protocol f(t) interpolates between adiabatic and instantaneous switching. For the 
Hubbard model this has been carried out by means of lowest order Keldysh perturbation theory [l8| . Due to the finite 
switching time, the dynamics becomes adiabatic for fast modes leading to /j, = 1 for fast modes (large momenta). 

For the ferromagnetic Kondo problem |14l [19[ a similar analysis of the crossover from adiabatic to instantaneous 
quenches has not yet been established. Since it posesses a nontrivial scaling flow, it is not clear what energy scale 
marks the transition from adiabatic to a sudden response of the electron system. In addition, there is no intuitive 
picture about the fate of the renormalized exchange couplings in presence of a time-dependent vertex. Applying our 
new time-dependent flow equation approach, we will show that the relevant scaling flow is just the instantaneous one 
(at least for the resonant infrared couplings which characterize the dynamics of the Kondo system at large times). 
This allows us to study the real time evolution for this non-sudden and non-adiabatic switching procedure within a 
non-perturbative scheme. Here non-perturbative means that this result cannot be obtained from summing a finite 
number of diagrams. The advantage of our renormalization scheme is to avoid the usual complicated reordering of 
the perturbation series labeled by the bare coupling J(t), by solving the truncated one loop RG-equations. 

Our motivation is therefore twofold: In the first part of this paper we introduce the new formalism and explain 
how it can be used for generic time-dependent many body problems. In the second part we apply this approach 
to study the quantum dynamics of the ferromagnetic Kondo model HJ subjected to a continuously switching on 
of the interaction. We calculate the mismatch factor [i from JS} as well as the long-time behavior of the impurity 
magnetization curve depending on the switching protocol. 



II. TIME-DEPENDENT FLOW EQUATIONS 

A convenient way to solve the time-dependent many-body problem would be the application of a unitary time- 
dependent transformation U (t) that makes the Hamiltonian 

H(t) = U(t) H(t) U\t) i U(t) (6) 

simple enough for solving the Schroedinger equation 

iBtMt)) uw =H(t)\1>(t)) uw (7) 

in the new basis \^{t))u(t) = U(t) \ip(t)). Now even for the simplest quantum models such as the driven two-level 
system, there is no general way to find a transformation which simplifies (J7|) considerably and therefore this approach 
is not commonly used. 

In this paper we will develop such a general approach based on a suitable extension of the time-independent flow 
equation method (method of infinitesimal unitary transformations) [l2l[l3j. Instead of performing the change of basis 
(|6|) in one step, it is convenient to use infinitesimal steps by solving the differential equation 



d B H{B,t) = [r](B,t),n(B,t)}+id t r]{B,t) 



(8) 
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with a suitable antihermitean generator t](B, t) and the initial condition posed at B = 0: T~L(B = 0, t) = H(t). Eq. © 
becomes equivalent with (jB]) for the following unitary transformation constructed as a -B-ordered exponential 

U{t) = U(B = oo, t) = T B fetf" «'»( fl .*)l . (9) 

This approach stays conceptually close to the time-independent formulation and becomes equivalent to it for time- 
independent Hamiltonians H. Now the key question is how to choose the infinitesimal generator r)(B,t) such that 
the series of Hamiltonians H(B,t) = U(B,t)(H(t) — idtjU' (B,t) eventually becomes diagonal during the flow from 
B = to B = oo. The suitable choice is to add an additional time derivative to Wegner's canonical generator for the 
time-independent case [131 | : 

r](B,t) = [H (B,t),H int (B,t)} - id t H int (B,t) , (10) 

Here Ho denotes the diagonal part of the flowing Hamiltonian and Hint its interaction part, H(B,t) — Ho(B,t) + 
Hint(B,t). Then one can easily show that the time-dependent interaction matrix elements satisfy an inhomogeneous 
diffusion equation with a negative source term 

(d B -d?)Tr[Hl lt (B,t)] =-2Tr[ V (B,t)rf(B,t)] <0. (11) 

This ensures that for B — > oo only resonant couplings Hf nt (B — !• oo,t) = % 3 -e _ *^" — remain in the interaction 
part of the Hamiltonian. Hence in general there remains a simpler time evolution problem in the B — oo basis with 
at most the complexity of the common rotating wave approximation. 

The main computational problem of our new approach is that solving the set of partial differential equations ([8]) 
and (jlOp is more difficult than the set of ordinary differential equations in the time-independent case. Notice that 
according to (|10p the parameter B has the dimension energy -2 , therefore our flow occurs in an energy-scale separated 
way by first decoupling matrix elements with large energy differences before looking at smaller energy differences. 
This will turn out to produce an RG-like scaling flow for time-dependent Hamiltonians. 

Next we illustrate this new approach for the weak coupling ferromagnetic Kondo model (JTJ), which will allow for 
a controlled truncation of the infinite set of flow equations in an expansion in the time-dependent running coupling 
Jkk'{B,t). The calculation proceeds along similar lines to the interaction quench in the ferromagnetic Kondo model, 
therefore we adopt the notation from Refs. [2[lj|- Details of the calculation will be published separately. 



III. TIME-DEPENDENT KONDO MODEL 

To leading order in the running coupling constants (equivalent to one loop order in the conventional scaling ap- 
proach) the flowing time-dependent Hamiltonian can be parametrized as 

H(B,t) =^£ fc 4 >CT c fe , <r + H ini ( J B,t) (12) 

with the interaction part 

H int (B, t)=J2 Jkk> (B, t) S ■ s k « (13) 

k,k' 

From ([8]) and (jlOl) one derives the key result for the time-dependent one loop scaling equation of the coupling constant: 

d B J kk ,(B,t) = -A 2 kk ,(t)J kk ,(B,t) (14) 

— ^2i{JkqAqk'{t)Jqk' ~ <J 'qk' Afeg(t) J \ q ) (n q — - J 
q V ' 

where we have introduced the operator 

Afcfc'W = (ek-£k>-idt) ■ (15) 

Here n q is the occupation number distribution of the noninteracting conduction band. In the sequel we will only 
consider the zero temperature case (the generalization to nonzero temperature is straightforward): n q = Q(—e q ) 
where we have taken the Fermi energy ep = 0. Notice that (IT4l) holds for arbitrary time-dependence and both 
antiferromagnetic and ferromagnetic couplings. The neglected normal ordered terms in the flow of the Hamiltonian 
are third order in the running couplings. 
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For time-independent couplings Eq. (fT4"f reduces to the standard one loop scaling equation of the flow equation 
formalism 12]. For the time-dependent case Eq. (fT4"|) has the structure of a complicated non-linear diffusion equation: 
An interaction that is initially localized in time spreads out during the f?-evolution. Notice that this does not lead to 
violation of causality since all the observables undergo the same unitary flow, therefore a time-dependent expectation 
value at time t does only depend on couplings J(t') for t' < t. 

Eq. (fT4"]) is best analyzed in Fourier space, 

Jfcfc'H = / dte wt J kk ,(t) . (16) 

The solution of the linearized form of (|14p gives the enveloping behavior of the flowing couplings which we later dress 
with the one loop scaling flow. The linearized solution is simple in Fourier space: 

J w {B,u) = e- BA "'^J(B = 0, W ), (17) 

with Afcfc'(w) = £fc — efe' — uj. The initial condition J{B = 0,w) is the Fourier transform of J(t) from (JT|). Therefore 
our formalism with the generator (|10[) eliminates the non-resonant Fourier modes first in an energy scale separated 
way in a direct generalization of conventional time-independent scaling theory. 

Now one can parametrize the full solution of ((TJJ very accurately (meaning: correct in one loop order) by the 
following approximation 

Jw(B,u) = e - BA "'^ £ _(B, W ) (18) 

where the flow of j e —{B, ui) is determined from the flow of Jkk' {B, u>) at the point A.kk' (w) = 0. The notation e^p- 

corresponds to the average energy of e k and ey: t-^r = [e k + £fc')/2. In (fT8|) the flow of the dimensionless coupling 
j e —.(B, uj) therefore effectively dresses the linearized solution with the one loop scaling flow. 
Substituting A/ c / c '(cj) = into (| 14|) yields the following scaling equation for j € (B,u>) 

dB ~ 2B P J 2tt [ } 

Xje+u/2(B, I + uj')j e _ u/2 (B, | - uj') 

Here the initial condition is posed at B = D~ 2 , where D is the band width of the conduction band, p is the density 
of states in the conduction band, which we assume to be constant. 

Specificially, for the time-independent case with Jkk'(B = 0,t) — Jo one finds j e (B = Q,ui) = 2ttS(uj) Jo. Since the 
flow at the Fermi surface determines the low-energy physics we use the identification j ep (B,uS) = 2it5(u)) J(B) and 
(1191) reduces to the usual one loop scaling equation for the infrared behavior 

^ = - P J 2 (20) 
dB 2B H y ' 

or in terms of the scaling parameter A = B^ 1 / 2 with dimension energy: 

dJ o 

d^K = - pJ ^ 
The solution shows the characteristic logarithmic scaling of the Kondo model, 

Jo 

1- ^ln(BD 2 ) 



J(B) = - ^ ■ (22) 



IV. CONTINUOUS SWITCHING OF THE EXCHANGE COUPLING 

In our model we are interested in a situation where the exchange interaction J(t) is switched on continuously over 
a timescale r. For convenience we choose the specific form 

Jw(B = Q,t) = f (l + erf (±-\\ (23) 
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with ferromagnetic Jo < 0: lim^-oo J(t) = and lim f _>. +00 J(t) = Jq. The Fourier transform of (|23p is 

2 2 

ie~ T u 

J kk ,(B = 0,u;) = Jo (24) 

UJ + jU" 1 " 

which permits simple analytic expressions in our subsequent calculations. For example the linearized solution (|17j) of 
(fl4l given the initial condition (|24l) takes the following form 

J kk ,{B,t) = ^a- B ^-*y) 2 (25) 



x 1 + erf 



-iB{e k -e k ,)+t/2 



T 



Notice however that the formalism can be carried through for an arbitrary switching protocol. 

Next we need to calculate how the one loop scaling flow (fl~9| dresses the linearized solution. Using the Dirac identity 
one can verify that the parametrization 

2 2 

ie~ T u 

j ep (B,u)=g(B,cj)—— (26) 
uj + 10 + 

in (fH)|) leads to the following scaling equation 

d -^ = ^ Pg (B,0) 9 (B^e-^/\ (27) 

with the initial condition g(B = D~ 2 ,uj) = Jq. One can verify that the shift tp ± uj/2 on the rhs of (fT9|) can be 
neglected in one loop order in the derivation of (|2"T)) . Eq. (|2"7| can be solved easily: For B < lu~ 2 the flow is to a very 
good approximation (meaning: correct in one loop order) just the static flow 

Jo_ 

1 _ £^l n (B£»2) 



9(B,u) = - pJo :,„^ (28) 



which stops once the energy scale B w uj is reached, 

g(B > oj- 2 ,uj) = ^ . (29) 

1- £^ln(£> 2 /cj 2 ) 

Remarkably the switching scale r does not enter in these equations. Notice that we find the expected weak coupling 
flow for the ferromagnetic case, which makes our calculation controlled by a small parameter. 
Summing up we have derived the following one loop flow of the coupling constants 

_ 2 2 

J kk ,{B,u) =g(B,u) l -^-r e- BA l^ (30) 
uj + iQ + 

with g{B,uS) according to (l2"5t and (|2"9")l for energies around the Fermi surface, (e^ + ey)/2 = ep. 



V. EVALUATING OBSERVABLES 



Next we want to evaluate observables making use of the diagonal structure of the Hamiltonian for B = oo. Specifi- 
cally, let us consider an observable O without explicit time dependence and an initial state |\&j(to)) given at time to. 
Then the expectation value of O at time t is 

(0(^^(^)10^)) (31) 

where ^i(t) follows from the solution of the time-dependent Schrodinger equation with respect to H(t). We can 
rewrite this in the Heisenberg picture 

(0(t))j = (*j(*o)|0(t, *o)|^j(*o)> (32) 

with 



0(t, t ) = U\t Q ) Pt( t) to) U(t) O U\t) P(t, t Q ) U(t ) 



(33) 
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and the time evolution P(t, to) evaluated in the simple B = oo basis: 

ft 

dt"H(t') 



P{tM)=% 



exp 



in 



(34) 



O(t,to) from (133)) is therefore the solution of the equations of motion for the operator O with the initial condition 
0{t = t ,t ) = O. 

The explicit evaluation of f|33[) proceeds similar to the forward-backward flow equation transformation scheme for 
time-independent Hamiltonians [2l|- First we derive 0(t) = U(t)OU^(t) by solving the differential flow equation 



Q(B,t) 
dB 



[r 1 (B,t),0(B,t)} 



(35) 



with the initial condition 0(B = 0, t) = O and the result 0{t) = 0(B = oo). Next we employ the time evolution 



0(Mo) d = PHt,to)d(t)P(t,t ) 



(36) 



which will usually be simple since % is diagonal. Then we undo the transformation and go back into the original 
basis, O(t,t ) — W(to) O(t,to) U(to) by solving the differential equation 



O(B,t,t ) 
dB 



[ V (B,t ),O(B,t,t )] 



(37) 



with the initial condition 0(B = oo,t, to) = 0(t, to) now posed at B = oo. The solution at B = gives the desired 
time evolved operator for (f3"2")l , 0(t, to) = 0(B = 0, t, to)- The main difference to the time-independent case [2l| is 
that here the forward and backward unitary transformations U(t) and U(to) will generically be different. The key 
advantage of our approach is like in Ref. plT | that there are no secular terms in time (i.e. errors that grow with t — to) 
since we can evaluate (|3"6")l exactly for diagonal Hamiltonians in the B — oo basis. 



VI. TIME EVOLUTION OF THE IMPURITY SPIN 

We want to use the above framework to derive the time evolution of the impurity spin expectation value starting 
from an initial product state with a fully polarized spin at to = — oo: 

= |FS) ® | t) (38) 

where |FS) denotes the zero t emp erature Fermi sea. For the unitary transformation of the observable O — S z we 
make the following ansatz [l2|, Il9l | 

S z (B,t) = h(B,t) S z + J2lkk>(B,t) (S+ S ^, k - SS+ k ) (39) 
fe.fe' 

where we neglect hig her order terms in order 0(J 3 ). The flow equations (1351) for the coefficients are structurally 
identical to Refs. [12|, [19J with an additional time dependence: 

d B lkk> (B, t) =h(B, t)A kk , (t)J kk , (B, t) (40) 

d B h{B,t) = J2~fk>k(B,t)A kk ,(t)J w (B,t)n k (l - n k >) 

k,k' 

It is possible to verify that one can use the solution (|30| for all values k, k' in these equations, the deviation only plays 
a role in higher orders. The solution for the Fourier transform at B — oo with the initial condition h(B = 0, t) = 1, 
Ikk 1 (t) = is readily obtained as 



2 2 

ier T w 



lkk ,(B = oo,u;)=g((e k -e k ,)- 2 ,0) ' (41) 

A fcfe / (w) (u) + iO) 

Time evolution (|3"B"|) just generates additional phase factors e^ £fe_£fc '^ t_t °^. Inserting the result back in (|3"?| and 
integrating back to B = yields 

7fcfc'(Mo) - g{(e k -e k ,)- 2 ,u = 0) (42) 



2tt A fefc , (w) (u + »0) 



'") 
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and 

h(t,t ) = 1 + \ 7 fefc / (t, t ) 7fc/ fc (t, t ) rik (1 - (43) 
fe.fe' 

plus terms that are third order in the running coupling. Eqs. (|42|l and (|43|l therefore constitute the solution for the 
operator equation of motion in the original basis and we can immediately read off the desired spin expectation value 

(S,(*))j= lim (9i\S t (t,to)\^i) = \h{t,t Q -> -oo) (44) 

Now one can easily work out 

e -i(efc-e k /)to 



Jkk'(t,t ) = g({e k -e k >) 2 ,w = 0) 



Cfc - £fc' 



x ^ e i(e fc -e v )* _ e-r 3 ^*-^/)^ (45) 
for t 3> max(r, -D -1 ), to <C — max(r, D^ 1 ) which yields 



2 4 ./o 

,-2r 2 e a 



X (1 + e~ 2r e - 2e~ r £ cos(ei)J (46) 

where J(e) = Jo/ (1 — pJo m (-D/ e )) is just the usual infrared running coupling constant (|22[) of the time-independent 
Kondo model. Eq. (|46|) is our key result which contains both the asymptotic magnetization for t — > oo and the 
approach to it as a function of the switching time scale r. It can be viewed as a systematic resummation of the short 
time perturbative expansion [2(| that incorporates both scaling effects and the switching protocol (f2"B"l) . Notice that 
while we have made the simplifying assumption t 3> r in deriving (|46[) , one can obtain the behavior on all time scales 
from (gl]) and (j43]l. 

Eq. (|46p shows the correct limiting behavior for adiabatic switching on (r — > oo) 

{S z (t = oo))/ = - - - / de = - + — + 0( J ) (47) 



2 4 7 e 2 4 

which reproduces the equilibrium magnetization for an infinitesimal magnetic field (22[. For instantaneous switching 
t = Eq. (|46|) reproduces the results from Refs. [HI 

{S z {t = oo)) I = l - + P ^ + 0{4) (48) 
The mismatch factor /i introduced in (J5)) interpolates between (l47|) and (|48|) as a function of the switching time r, 

<&(< = oo))/ = i + /i(r) ^ + 0( J 2 ) (49) 

and can be derived from (|46"1) 



(50) 



It is depicted in Fig. Q] and shows the expected crossover from 2 (instantaneous quench) to 1 (adiabatic quench 
equivalent to equilibrium). Notice that /i(r) starts to deviate from 2 when r becomes of order the inverse bandwidth 
since then there are fast modes with e > r _1 for which the quench is no longer instantaneous. However, the approach 
to the equilibrium value 1 is very slow, 

^ >>D ^ = l+ l- P jl H rD) (51) 



For example the value jii(r) = 3/2 is only reached at the nonperturbative time scale t = D 1 cxp [ — j^jjjj 

Likewise the approach of the spin expectation value to its asymptotic value is also logarithmic and (remarkably) 
independent of r as can be derived from (|46j) 



for t 3> max(r, D 1 ) 



W' = ^ = »»'-5liiPl (52) 
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FIG. 1: Mismatch factor fi(r) from (O for pj = -0.1. 



VII. CONCLUSION 

We illustrated how a generalization of the flow equation method can be used to study time-dependent Hamiltonians 
with a nontrivial scaling behavior. In particular for the time-dependent ferromagnetic Kondo model we calculated 
the mismatch factor /i which measures the "non-adiabacity" of the quench. Different from previously studied models 
like the quench in the Hubbard model (l5l . [l6l [l8| here the approach of \i to the equilibrium value is only logarithmic 
and one needs to quench very slowly to reach the adiabatic limit. 
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